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ABSTRACT 


The propagation of intense radiation in a magnetized 
plasma is investigated theoretically. The presence of 
intense radiation in a plasma can alter the plasma density 
profile in the region of the beam. This alteration can be 
caused by the ponderomotive force or by the heating of the 
plasma by the radiation. We have calculated the equilibriun 
density profiles and hence the dielectric € , for each 
mechanism. The wave equation was then solved using a WKB- 
method and differential equations were derived that govern 
the propagation of the beam in the plasma. It is shown that 
the behaviour of the bean propagation is determined by a 
competition between diffraction of the beam and the focusing 
tendency of the beam produced by the altered density 
profile. The nature of the solution. to *he wave equation is 
discussed along with applications to the laser solenoid 


fusion concept. 
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1. INTRODUCTION 


1.1 THE LASER SOLENOID APPROACH TO FUSION 


Studies of heating of magnetically confined plasmas to 
thermonuclear temperatures have led to the investigation of 
several possible heating schemes, one of which is the laser 
- solenoid approach proposed by Dawson et al (1971). fTfhis 
scheme consists of a long slender plasma column confined by 
a solenoidal magnetic field (100 - 500 kG) and intense, long 
wavelength laser radiation (ie.10.6m CO, laser) directed 
axially so as to heat the plasma colunn to thermonuclear 
temperatures (™~108°%K). This simple geometry has the 
advantage of being relatively stable and the heating source 
(CO, laser) holds considerable promise in being scaled up to 
the reguired output. However, the physics of the laser 
plasma interaction must be carefully analysed both 
theoretically and experimentaily before large scale 
solenoids can be constructed. In this thesis, we consider 
theoretically one aspect of this interaction; the bean 


trapping or seif - focusing problem. 


1.2 REVIEW OF PREVIOUS WORK ON BEAM TRAPPING AND SELF - 


FOCUSING 


Steinhauer and Ahlstrom (1971) showed that a laser 
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propagating in a theta pinch would be refracted out of the 
plasma (ie. become untrapped) after propagating a distance 
of a few beam widths because of the density maximum on axis 
produced by the pinch. As light tends to bend in the 
direction of increasing refractive index N, and since N 
depends on the plasma density n, by N=(1-n/n™) (where n® is 
the critical density for the laser radiation), then if the 
density had a minimum on axis, the laser light would tend to 
bend (refract) towards the axis of the solenoid instead of 
refracting out of the plasma. Using a parabolic density 
profile, Steinhauer and Ahlstrom (1971) employed ray optics 
to show that the rays became trapped in the plasma and 
oscillated about the axis as they propagated along the 


plasma column. 


Humphries (1974) analytically solved for the waveguide 
type modes in a plasma with a parabolic density profile and 
determined the width of the smallest fundamental mode. Mani 
et al (1975) solved the wave equation with a parabolic 
density profile and a Gaussian intensity profile and 


produced solutions which showed the beam alternately 


focusing and defocusing in agreement with ray optics. 


In all the above analyses, the plasma was assumed to 
possess an arbitrary density profile, however it is known 
that the presence of an intense laser beam can produce a 
density minimum on axis either by the ponderomotive force 


(Hora (1971)) or by the heating of the plasma (see Burnett 
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and Offenberger (1974) for example). Thus the effect of the 
iaser on the plasma can in turn affect the laser beam itself 
(by causing focusing) so that a self - consistent treatment 
is necessary. Self - focusing is the result of the laser 
producing its own favourable density profile either by the 
ponderomotive force or by heating. TOC SRE Peat. Se. aE oc 
consistent treatment of the ponderomotive mechanism in an 
unmagnetized plasma was given by Sodha et al (1974) and Max 
(1976) where oscillatory solutions were found. The heating 
mechanism was investigated analytically by Sodha et al 
(1974b) for an unmagnetized piasma and numerically by Feit 
and Fieck (1976) for a magnetized plasma. Both of these 
investigations showed the alternate focusing - defocusing 


behaviour of the beam propagation. 


Several experiments have observed beam trapping and 
they are listed in the references. However, they only 
demonstrate this effect over short distances and for 


temperatures much less than the thermonuclear regime. 


1e3 NATURE OF THE PRESENT WORK 


In this thesis we extend the results of Sodha (1974b) 
and Max (1976) by including a magnetic field inside the 
plasma. The procedure used is to compute the steady state 
dielectric of the magnetized plasma (Chapter 2) and then 


solve the wave equation (Chapter 3) using the WKB method of 
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Akmanhov et al (1966) and obtain differential eguations 
which describe the focusing and defocusing of the laser bean 
in the plasma. We then investigate the nature of the 
solutions to these equations in Chapter 4 and discuss these 
results and their ramifications for the laser - solenoid 


fusion scheme in Chapter 5. 


Wwe consider in this thesis, two mechanisms that will 
produce the favourable density profile for focusing. These 
are; the ponderomotive mechanism and the heating mechanisn. 
In section 2.1 we compute the steady state dielectric € due 
to the ponderomotive force. Since ion motion is involved 
(ie. plasma is pushed radially outwards from the laser 
bean), pulse times are required to be larger than the 
acoustic transit tine Ta=%, in order to establish the 
egilibriun (where c is the ion sound velocity anda is the 


1fe beam width). 


In section 2.2 we extend the work of Sodha et al 
(1974b) and use the kinetic formalism of Sharkrovsky et ai 
(1966) to compute the dielectric € for the case where 
heating causes a redistribution of the plasma. As the 
results of this theory turn out to be quite restricted, 
especially to low intensities, we will denote these results 
as "weak heating" resuits. The weak heating results are 
extended in the next section so as to include the effects of 
ion heating and thermal conduction where we use a simplified 


fluid theory approach to compute € These results will be 
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denoted as "strong heating" results to differentiate then 

from the results of section 2.2. In both the heating cases, 

the heating time, Te3nT (woere Tf; is an ev, 4, is the 
dy 

absorption coefficient and I is the laser intensity in 

Watts/cn*) must be larger than the acoustic transit time in 


order to achieve the desired equilibrium. 


To simplify the analysis we treat the laser radiation 
as being linearly polarized which is valid for w >> 
(where @ is the laser frequency and ff, is the electron 
cyclotron frequency), as well as consider the radiation to 
be monochromatic. We also neglect the axial plasma dynamics 
and thereby consider only the radial equilibria, which means 


that our treatment is not fully self - consistent. 


Cgs units are used throughout except we adhere to the 
convention of expressing temperature in electron volts (eV) 


and power in Watts. 
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2e SELF-CONSISTENT CALCULATION OF THE NON-LINEAR DIELECTRIC 


IN A MAGNETIZED PLASMA 


In this chapter we compute the dielectric € of a 
magnetized plasma in the presence of strong electromagnetic 
radiation. We begin in section 2.1 by employing a fluid 
model for the case of a collisionless plasma where the 
ponderomotive force is responsible for altering the plasma 
density profile. In the next two sections we use kinetic 
and fluid theory respectively to calculate an expression for 
the dielectric where heating of the plasma by the radiation 


field alters the density profile. 


2-1 COLLISIONLESS PLASMA - PONDEROMOTIVE FORCE 


We begin by treating the electrons as a fluid and solve 
for the motion of the electron fluid in the presence of 
strong electromagnetic radiation (ie. laser). The equation 


of motion cf the electrons is 


ode _ —y 
{2.1} yD+(B-V)V= 2 €- 2 UxB - V Pe 
m me Cee 


where P, =“NeTe- We solve {2.1} by separating oscillatory and 


non-oscillatory quanties; 
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where VU, is the unperturbed fluid velocity and VU is the time 


varying perturbation induced by the radiation field. The 
total electric field € is composed of the time varying laser 
field EF and an electrostatic field ‘F,,due to any charge 
separation between electrons and ions. It is assumed that 
no external d.c. electric field is applied. The total 
magnetic field B is composed of the laser field contribution 
B and a constant, externaily applied d.c. field, B=(0,0,8). 
All oscillatory quantities are taken to vary as fot where 
W@W is the frequency of @ the’ Pradiation™ field, * so’. that 
substitution of {2.2} into {2.1} yields the zeroth order 


velocity (for as 2>>Q2=e2B2/m2 ce, and E>>T/feL where Lisa 


typical scale length), 
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We restrict discussion in this section to the collisionless 
time scale, so that if no collisionless heating takes place, 


we have 
ae — 
{2.4} VPe mS ays 


Using {2.3} and {2.4} together with the Maxwell equation 
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in {2.1}, we obtain after averaging over fast time scales 
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The second term in {2.6} represents the non-linear 
(ponderomotive) force due to the oscillatory field. tes 
to be understood that the density n, in the above equation 


is the time averaged density <n,> . Using Ampere's law 
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{2.6} becomes, 


AGE Vne : efe,- € (F-V)E‘+ Ex¥x E*] (Vx? DP, 
Ne ale: 2m 0 Te 477 


ae => aw ow =— rte — ee 
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{2.10} (Vx BL)x Roe GRE, es 

in {2.8}, we obtain 
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Ne £ 4h w* na Te SH 


where we have assumed that the external magnetic field B, to 
vary only in ‘the “radra?}l’ “dinectrvon. Equation {2.11} 
represents the equilibrium equation satisfied by the 
electron fluid in the presence of a radiation field. The 
second term on the right-hand-side (the ponderomotive force) 
Can cause a redistribution of electrons thereby producing 


radial density and magnetic field gradients. The ions are 


Mies oF Ris 
Oh a) ve 
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coupled to the electron motion through the charge separation 


field £,,. 


If the ponderomotive force and the VB2 force have scale 
lengths greater than a Debye length, then we can neglect the 


electrostatic field E over these scale lengths and write 


{2.11} as 
V VE’ VR 
{2.12} n, -~= Oo Vimeo! 
5 { og ial 


where = €) a BS 
e Fon, to e 


We solve {2.12} in the infinite conductivity linit, 


ie. assuming frozen magnetic field lines, such that 


{2.13} le B 
B 


8 


where 2,,B, are the density and magnetic field far away from 


the laser beam. The soiution of {2.12}, using {2.13}, that 


satisfies the boundary condition; 1—?n, as E=- 0 is 


ae 
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{2.14} dni. ta Zo (Het) Srete F 


where 
y= ne /Nn, 
{2.15} (2.= tn, Te 
Rie 


For sufficiently high magnetic fields, BX<<16 and {2.14} 


becomes 


{2.16} ne die x, BEY 


Thus in this limit, the equilibrium density profile is 
parabolic in E as compared to the exponential dependence 
(Max (1976)) found in the field free case (B=0). The axial 
magnetic field can therefore reduce the plasma depletion due 


to the ponderomotive force, as expected. 


With the knowledge of the equilibrium density profile, 
mn: we can write down the steady state dielectric valid for 


uf >>92, and collisionless time scales as 
€ 
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{2.17} Ge Sn ae &(€) 


ees 
r z 
where Wp = ey wt nem, 
inter a 
{2.18} Cpa) 


PCE)= Hp ( 1 — m(E)) 
re he 


The non-linear term, @(E2), of the dielectric vanishes for 
E2 —-? 0 as the plasma becomes a linear mediun Cer76) with 
no strong electromagnetic fieid present. For the case of 


the ponderomotive force in a strong magnetic field (8,<<1) 


$ becomes 
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This expression will be used later in the solution of the 


wave equation and the subsequent investigation of self- 
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focusing. 


2-2 COLLISIONAL HEATING - KINETIC THEORY 


In this section we compute the steady-state intensity 
dependent density and temperature profiles where collisional 
(inverse - bremsstrahlung) heating of the electrons is taken 
Tuto. account. As this process involves electron-ion 
collisions, the time scale required to establish an 
equilibrium is much longer than the acoustic transit time 
Cae so we are dealing with a longer time scale process 
than in the case of the ponderomotive mechanisn. We 
restrict the analysis to intensities such that the 
ponderomotive force discussed in the previous section may be 


neglected. 


We begin by solving the Boltzmann equation; 


{2.20} a, 5 + V+ d29 a OF (B,xv)- 2 52 IF/f¢ 


> end _ ie - : po. 
where a=-eE/m,, 2=cB/m and Jd f£/St is a’ suitable collision 


term. 


The solution of {2.20} wili follow the formalism of 
Sharkrofsky et al (1966) whereby we expand the distribution 


function , f, in terms of spherical harmonics; Yeme) ie. 
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i) = E SoM 


ms 


= $ 


+ aye. TOBY =) +5 sin ae + 5, sin 6 Sengt.-. 


pet} Sle eae Ue: 


where f,=f, (vw) is the symmetric portion of the distribution 
function and 
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is the: first order, asymmetric correction to the 
dzstrapuciron ~Luncc1 on. Expansion {2.21} was used by Sodha 
et al (1974) in the investigation of self-focusing in a 


field-free (8=0) plasma. 


According to Sharkrofsky et al, expansion {2.21} is 
useful if; 
(G) in one wave period, the change in the electric field 
seen by the particle is. small. This condition can be 


represented by 
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2 
{2.228} Ve\ cc | 
ALS 


where v, is the electron thermal velocity. 


(ii) the electric field be such that the electron quiver 


velocity, V=eE/aw, is smail compared to UY yr Lee 


ae a i 
{2.22b} ey Ze | 


Substitution of the spherical harmonic expansion {2.21} 


into the Boltzmann equation {2.20} yields; 


1. The 1=0 scalar (density) equation 


= —_ 
2.23 10.-5,4+.3 (wa-§ )- ; a = 
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where Pis the electron-ion collision frequency. 


2. The i=1 vector (mometum) equation 
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o> 


= es mete 
2.20) 2,54 vdp5,+ dH + Rexd, + VF 


Due to the properties of the spherical harmonics, 
ensemble averages of scalar and vector quantities are 


defined as 


YePanse <Hys te PE v'd 


{2.25} 
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Q=ayB > <A> 47 (@ : 


In the presence of time varying fields, we assume a 


‘wt ; 
time dependence of ade and expand the distribution 


function as a Fourier series; 
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From {2.25}, we see that the density is given by 
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Ckut ear, 
2204 n= 47 > e > vidv 
K Oo 


where the real part is to be taken. We can define 


DSi as6 aS 


} ‘kot; k 
k t kwt t s be 
{2.28} n € = 47 € 5 valv , k=0,1,2, 
6 
where the igs Kol 1alevassoci1a Cede) with | a.c. Space charge 
effects. These effects can be neglected for the case where 


thermal velocities are much less than the phase velocity of 


the e.m. wave (U,<<b/k~c) and we thereby set 
iG 
3 
{2.29} N= n= on{ 5 v°dv 
° 


and n* =0 EOD) Koel. From {2.25} we see that the average 


velocity is 
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vw 
ckwt ( >k 3 
{2.30} p= ane fa5 Meee 
3n - 


In the following treatment, we simplify the analysis by 
considering only the d.c. and fundamental enone terms 
(k=0,1) in the expansion for f. This is justified when the 
highest frequency of interest is the oscillation frequency 
of the e.m. wave since higher harmonic terms will have 
frequencies considerably removed from any other frequency of 
interest. This impies the frequency ordering by 2 >>y? RP oly? 


etc. 


We substitute the Fourier expansion {2.26} into the 
Boltzmann equation {2.20}, retaining only the k=0,1 terms 
and obtain the resulting density and momentum equations, 
which will each split into oscillatory and non-oscillatory 


parts. The 1=0 (density) equation yields; 
Vd 
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{2.31b} Py _v 


The momentum equation yields; 


o .3° ry = | 1 
{2.32a} Vidz t+ 4 On, Fo +t he a At oe ) + 


m4 = ie Oh 
eet 3, I= © 


x . a O i] 
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Daas A ‘ Po = 


In the above equations, “a9=-6E9/n, where E9 is an 
electrostatic field (ie. due to charge separation) and “al=- 
eE1/n, where E1 is the amplitude of the e.m. wave. The 
neglect of a.c. Space charge effects implies that n1=0 
which, from {2.28} implies that f= We therefore set £.=0 
in {2.31} and {2.32} and solve them for f°,f,°,and F', take 
the external magnetic field along the z-axis only, then from 


{2.32a} it can be shown that 
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| | 
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SUbStatuciony Of {Ze s3})e.and {2e35)) nto. f2.32a}, wi: 
yield an equation for the symmetric part of the distribution 
function a5 (v)- We write the gradient operator as V =2,+d, 
_where the directions are with respect to the external 
magnetic field ge If the following assumptions are made, 

1. linearly polarized e.m. wave 

2. E°<X<E! (or scale length >> Debye length) 


3. W2 >> 2,2 >>y2 


the equation for ia becomes 
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In order to Simplify the above equation, we 


dimensionaliy investigate 


the ordering of terms in {2.37}. 


Taking the scaie length for perpendicular gradients to be L, 
we find that 


Uv 
term (1): mu ~ (4) we ar 
3 L* A] 3 
term (2): G5 poe elt sx) 2 ie 
ce ay ae Z}] 3 


* om . ° o A Son 
term (3): Vi =o(assuming f> uniform along B,) 


“a2 


- 


ve — evodé § st piziguhe og Bren + 
+40, 83 ai sete? Le jaidhilion Pre es nplaaeee , eiiz 


yi od 028 a#ackSony tee i a le ton atpost siece neg bal 
Dy 
Ow 


23 


term (4): by symmetry the cross product vanishes 


o 
fern (5): sie Ted ep 
3 
where Jame and ae is the electron Larmor radius. Thus 
Py thre 


the ordering of the surviving terms is 


| ae Rp (2 ree 
L (ak 


Therefore the neglect of terms (1) and (2) in {2.37} 


requires 


(Ss \ << ; 


Since {~ 20e/m,", we can write this condition as 


2s | 
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Thus the ion Larmor radius must remain much smaller than any 
scale lengths of interest, and this implies that radial 
thermal conduction losses from regions of size ~ L will be 


negligible. 


Assuming that {2.39} holds, the equation for fy. is 
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The first integral of {2.40} is 


{2.41} Es (2w Daa coe Ve g 
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Since fand Det —> 0 asu-—?cO (more rapidly than ur 


—pPoo) then as {2.41} must hoid for all u, C must be zero to 


satisfy the limit u—»*e. Thus we must solve 
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that 
3A 
A= n ( Me 
aT Te 

so that 

ry Ww 2 
{2.43} aE =i We / ine exp (eee) 

277; 2Tz 


The above distribution function is Maxwellian with an 
effective electron temperature Tp =T,(1+%, E?) . Using {2.43}, 


ran and £,| may be computed. 


We now derive an equation for the electron density fron 
the above results considering only the steady-state, which 
is defined to be when the zero-order (d.c.) current, Jo. is 


zero (ie. y°=0, see {2.30}). Setting 
—* g 3c? 
{2.44} AE ne) e | ay 5 A hat imo: 


and using {2.43} and {2.33} in {2.44} yields (see Appendix 
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{2.45} aiGay U3) S = 2n.eF 


This is the equilibrium equation for electrons. 


We can repeat the entire procedure for ions and solve 
for eye etc. and obtain exactly the same equation as {2.45} 
except that Tp, would be replaced by T;=T, (a constant), 
because of the assumption that the ions are not heated by 
the. e.m. field. Thus /for ions) the equilibrium equation 


would be (e—9 -e) 


{2-46} T, oy n= 2n,e i 


We solve for n, by assuming that scale lengths are much 
larger than .a Debye length so that we have the quasi- 
neutrality condition XD: , and eliminate E° from {2.45} and 


{2.46} to obtain 
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and 


{2.48} meetin = (VC eNLaVA De 

Wo 
is the equilibrium density profile in the presence of the en 
field in a magnetized plasma. In the unmagnetized case, 
collisions play a dominant role and Sodha et al (1974) have 
found the equilibrium profile to be, by a similar kinetic 
theory procedure, 


- S/ 
{2.49} = (4% EY ) - 


Thus there is a significant increase in the depth of the 
density well in the B=0 case as expected. It must be 
stressed however, that {2.48} is not valid for arbitrary 
magnetic fields (B,) and intensities (E@). The assumptions 
used to derive {2.48} restrict the values of these 


parameters and we review these assumptions below. 


The conditions for validity of {2.48}, along with 


assumptions about symmetry, are 
1. wee tUt/L) 2 << 1 
2: Ree ae << 1 


3. p2<<92<<ey2 
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Conditions 1 and 2 are required for the validity of the 
truncation of the expansion of the distribution function 
whereas the last two conditions simplified the analysis. 
Conditions 2 and 4 are the most stringent and table 1 shows 
the values of parameters required for validity for various 
initial temperatures. We assume CO, radiation, and the 


scale length L to be the beam size, I~0.1 cm typically. 


The results obtained here (density and temperature 
profiles) have been compared with those from ae one- 
dimensional MHD code (Milroy (1976)) and found to agree weli 


(see Figure 1). 


We now compute the non-linear dielectric explicitly 
using the formalism of the expansion scheme of Sharkrofsky 


et al (1966). From the definition of the a.c. conductivity 


One we have that 


' ' eo 30! 
{2-50} a =n E. = tre vi 5, dv 
0 


a 
Substituting for £4 from {2.35} we see that G.1! can be 


written as 
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Putting the above expressions to use in {2.51} together with 
the expression derived earlier for ix gives, after 
integration 


= Wt (p Send 7 m Ff 
{2.52} Pee =i oe nO = Ge Pe FP FeRI +e, £7) 


Bonk J SVR A, Cs 


|! « rt y ns 
= | 
s . tig. ' 2a? os a: sa0hl, whee - 
- i | 


= Ps tet 
fz 


\ ! : 
(2265300 jo¢ Shr Sp ye Sauee See Pee 1 te NeLe 
BUF tm (14%,E)™ am coe 


and ali other components vanish. Now the components of 


dielectric tensor are 


, / 
{2.54} Ey, = ee eid aa ox 


and this becomes, “using {2.52} and {2.53}; 


é, Cee 
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o) o & 
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This non-linear dielectric will be used later in the 


solution of the wave equation. 


2-3 COLLISIONAL HEATING - FLUID THEORY 


In this section we try an approach that is less 
restricted than the more rigorous kinetic approach used in 
the previous section. For example, the the kinetic approach 
dra not ‘allow for ion ‘heating or for radial. thernal 
conduction whereas the following theory removes these 


Eestrictious: 


We start from the static fluid equations where inertial 
terms are neglected (ie. 2 =0) which is a valid assumption 
if we are solving for the equilibrium state (ie. for long 
time scales). The momentum equation becomes a pressure 


balance equation (B, is the external magnetic field), 


| 2s 25 
{2.59} ni (tea Ve) te Bees 8. 
va 57 


The electrical conductivity of the plasma is assumed to 
be infinite, valid at high temperatures, which gives the 


frozen field condition 
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{2.60} y = ne é 
No 


° 


The electron and ion energy equations are; 


fot yes (Le Rs ¢ 2 end. e's lem ac 
StF ee 3" 


2 
{apo}, = 
3 


In the above equations, I, is the laser intensity (W/cm2), 
his the inverse-bremsstrahlung absorption coefficient, Ree 
are the electron, ion thermal conductivities and @% is the 


e-i equilibration time. 


We accomodate the ion heating by assuming a fixed ion- 


electron temperature ratio for all spatial points, ie. 


ata lies Cs Oo <es| 


This assumption is related to the temporal structure of the 


laser profile (Vagner et al (1975)). 


Using {2.60} and {2.63} in the pressure balance 


equation, we obtain an expression for the density, 
t/ 
eS 2 L 
Ske Ke 140) P-/4 + | (146) 8.7 


where we have defined 


{2.65} B.= Yan. le 
Eo 


For strong B-fields, OX<1 and {2.64} becomes 
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{2.66} OOF.b Fa (140) Po JA_ 


We now have % as a function of electron temperature, so 
we must use the energy conservation equations to solve for 
T, in terms of the input laser power. AQGANG SP HZ2e6 1g. eto 


{2.62}, using {2.60}, we obtain an equation governing toe 


{2.67} 2, a 4 nannies + Op Ie OD, Xe’ $ T, ky aa. 
he 6X,° 


where we have assumed that the electron radial thermal 
conductivity is much smaller than the ion thermal 
conductivity in magnetic fields of interest. We can rewrite 


“4 
therecaX:. term in {2.67} by noting that 2.=C/T, , therefore 


OreK =- CO MTe = =-L2%0% and {2.67} becomes 
K 20st ee 
{2.68} 


The laser intensity is taken to be Gaussian with bean 


width % so that we have 


{2.69} Ty: es ky exp (- 7 ono) 


where, I. (w/cn*) 7 is the on-axis laser “1atensity. The 
inverse-bremsstrahlung absorption coefficient and ion 
thermals conductivity si(Milcoy (1976) ?Gandwa Spitzer: 44(1967) 


respectively) are 
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wie = 
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Ip Fz 
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{2.71} Men Poe e (an /Ae 7, ) ta 
Lv 


where r\ is the laser wavelength in microns and 1InA is’ the 


Coulomb logarithm. 


Using {2.69} to {2.71} the last term in {2.68} can be 
written as 


{2.72} Dehy eee) Aaa (= r?*/2 022) 
| o 2, 7 


where Z=1, and 


=F eat/ 
{2.73} Awy= £. 4x10! A (ps) R mir) © Leag 


We solve {2.68} by assuming the form of the electron 


temperature profile to be Gaussian, i.e. , 
x 2 
{2.74} T. = tetera ae, )+ Te 


This assumption is reasonable since the laser intensity is 
Gaussian and a form similar to {2.74} was obtained in 
section 2.2 using kinetic theory. Equation {2.74} has three 
parameters to solve for; the width 9%, the maximum on-axis 


temperature T,, and the temperature far away from the axis 
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ins To simplify the analysis, we restrict ourselves to the 


region near the axis and assume T, >> Tf 


> «€©SOd«hthhat we may = use 


the approximate form of the temperature profile 
2 ru 
{2.75} oe a exp (- r A25- ) 


Substituting {2.75} into {2.68} yields the equation 


76) Pe yt TACO opp (1/2 oy) — 2 = (© 
2 


e 


We set r=0 in {2.76} and obtain an expression for the on- 


axis temperature Tye 


2 


UV 
{2.77} Te > a) Alo) 
2 


Differentiating {2.76} twice with respect to r and then 
setting r=0 yields a relation for the temperature profile 


width -; {see Appendix B) 
{2.78} Co, 


Thus the temperature profile is wider than the laser 
intensity profile by + 5/2, which is a more realistic 
situation than that found in section 2.2 where the widths 


were the same. Equation {2.78} shows the effect of a non- 
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zero radial thermal conductivity. Note that {2.77} and 


{2.78} are only valid for T, >> To. 


Using) {2.77} and {(2.78}/an) {2.64} we ’canm ‘write the 


density profile as 
2 WA, 
a.79) Beli earzto BexpC rr) —s1bar Cone) 


~ 14 ~I/ 
where a S. 23X10 & on. ri Cite) © 7 
Bs 


Setting r=0 in {2.79}, we obtain an expression for the 
on-axis density ratio (0), 
{2.80} 0) = ( P+ an, Ess) 
Equation {2.79} can now be written as 


Va 
{2.81} Ail) = [1 + Xir)] - Xiry 


where 


{2.82} Xery= MSE. exp (-0'/ee;*) 


See ee eee 


/ 
(iFiarE i 


The real part of the dielectric can now be written as, 
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where 
ru 
(2.84) Pexy= Wp (1- Bx) ) 
Cs 


The on-axis (electron) temperature f, and density ratio 


eb are plotted in Figure 2 for various laser powers. 
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3. SOLUTION OF THE WAVE EQUATION WITH A NONLINEAR 


DIELECTRIC CONSTANT 


In this chapter we solve the wave equation uSing a WKB 
- solution which was developed by Akmanov et al (1966) for 
the study of e.m. wave propagation in liquids and crystals 
with nonlinear dielectrics. As waS mentioned in the 
introduction, we assume the frequency ordering W2 >> 2,2, 
such that we can neglect the non-vanishing off-diagonal 
elements of the dielectric tensor because these terms are of 
0(2,2/w 7?) and hence very small. We thereby study linearly 
polarized solutions of the wave equation with a dielectric 


of the forn 


{3-1} e=-€,+ $ +e €’ 

where 

323) g = Wee (1-%) 
wt 


and the effect of the magnetic field is to alter the density 
profile. The imaginary part of the dielectric is assumed to 


be much smaller than the real part since W?>> ~pé. 


The wave equation is 
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Equation {3.3} can be analyzed by employing the concept of 
Slowly varying amplitudes (Akmanov et al (1966)) whereby a 


} JON ergi ag (or di verging) beam can be represented by 
- 
- e E 


where € is a unit polarization vector and s is a parameter 
which accounts for the difference between the actual bean 
and a plane wave state. The difference is due of course to 
the action of the nonlinearity or to diffraction. 


Substitution of {3.4} into {3.3} yields, 


zu 


° zu im 
{3.5} Cv, + eet ae kd, - ko + wie )Ecue,s)= Oo 


_ irs 2° 
where we have negiected V[é- ¥¢] compared [OK € which is 
valid for W2 >> Wp 2. We solve {2.5} with cylindrical 


symmetry (V,2=0,, +r-19) and by neglecting dag compared to 
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k2f. This latter assumption implies that the scale length 
of variation of the amplitude of the electric field along 
the direction of propagation must be greater than the 
wavelength of the radiation field (i.e. Slowly varying 
amplitude assumption). Foliowing Akmanov et al we take the 


WKB-solution for the amplitude to be 


~t k S042) 
{3.6} as Ecr2zye 


where s=S(r,z) is the addition to the eikonal which accounts 
for the nonlinearity in the dielectric and E(r,z) is the 
amplitude of the e.nm. wave. Substitution of {3.6} into 
{3-5}, with a,,£ << k2€ yields an equation that can be split 


into real and imaginary parts as 
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Equations {3.7a} and {3.7b} are similar to those obtained in 
the study of lens-like (linear) media whose solutions are 
cylindrical or spherical converging (or diverging) waves. 
The difficulty in the above equations is that the amplitude 
E(r, 2) and the eikonal StE,z)™" arevscoupled*==by* the 
nonlinearity. Following Akmanov et al, the solution of 
toe Fo Loresstr, 2) Wier petegencraiized==torcy Lindcredal Yor 


spherical waves with variable radius of curvature whereby we 


take 
{3 . 8} £Cv 2) = irs Pts 7 Ce) 
2 


WneLve ()(Z) 915 i1nterpereted as the radius Or curvature, “and 
¢ (2) is the change in the phase of the eikonal due to the 
nonlinear medium. The solution for the amplitude will be 


obtained by employing a Gaussian- shaped ansatz, 


r z u = g 
ey ENS BS exp(-eYat stn) Fee 
5 °(2) 
where a is the 1/fe beam size at z=0, £(z) is a dimensionless 
beam-size parameter, F(z) accounts for the absorption by the 
medium, and E2 is the on-axis field intensity. Equation 
{3.9} assumes that the laser intensity profile remains 
Gaussian at all times, however the on-axis intensity 
E2/f2(z) and bean width af(z), change as the beam propagates 


through the mediun. 


In the following dicussions we will neglect the effect 
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of absorption over length scales of interest. This will be 
valid if the plasma is sufficiently hot so that classical 
inverse-bremsstrahklung absorption is weak, i.e. the plasma 
aS /Opticallysthin ,over, Jengthsa scales, of interest. The 
neglect of absorption implies that the imaginary part of the 
dielectric (i.e. €° in {3.1}) can be neglected compared to 
the real part. The discussion is exact for the case of the 
ponderomotive nonlinearity if no collisionless absorption 


process occurs. 


We substitute {3.8} and {3.9} into {3.7a} and {3.7b} 


Watour (z)=1 (1.e..) no absorption) and obtain” from {3.7b}, 
-| 
{3.10} Ble>= § >, 5c2) 


and from {3.7a}, using {3.10}, 


As {3.611} must hold for all r and Z, in particular r=0, we 


set r=0 in {3.11} to get an equation for the phase @(z)e 


ru 
fst 2t e5 ((z)= -/ 4 + Me ge ‘é 
Ws a cope 


where 


A At oy ie, 
pik sulaeayn aan 


oy bos dqneg. snbehis lala, in 


od Yo baa odd 203, MaKe a) co 
noi sgzonda sn i ce en 


4 


fat if} Bap is¥.t} 
vot ai a aissdo ae (folses n 


ee a cng <6 


; 7 - val a’) c Dei : 


. - | “Yer calek « LaF Ey Sie “Babole 
| | yo: ipa lige "—¢9 = 


‘2 “(pO M62") Mie Hod aia 
Thilaka a = 


eh os gsivolszey at 2 Soe. ils 308 blow Sas. cc. 
: Pesan Jebrsbs “aed set aouships as top of ur wa > 


2 & i 2 


; anne 


43 
: py rv zu 
3013} [T= c?kK —- UW + Wp 


The quantity F represents the nonlinear wave number shift 


due to the change in the on axis density. 


In order to obtain an equation for the beam width 
parameter £(z), we differentiate {3.11} twice with respect 
fomlre sande iseth-r-0,. andmeuse 13. 127 tosnget,sarter come 


manipulation, 


{3.14} kas a oF F 
9,2 3¢2)= 
/ 
Chee ea bs SS Hoye Yee) 
{3-15} co nt os 


where {3.14} is for the ponderomotive and weak heating cases 
discussed in 2.1 and 2.2 where ri] was a function of E. 
Equation {3.15} applies to the strong heating case discussed 
in 2.3 where $=2(%. In the above equations, the prime 


denotes differentiation of r3 With respect to its arguement. 


Equation {3.14} and {3.15} describe the propagation of 
the laser beam through the plasma over distances where 
absorption is not important. We will refer to the 


differential equation for f(z) as the self-focusing equation 
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for the particular mechanism involved (i.e. ene eoueE Ves 
weak or strong heating). The beam will converge (i.e. the 
effective beam-width af(z) decreases) only if the slope of 
PAZ) 25 a Gecreasing Lunctrvon ot “2, i.e: 2, £ (z) <0. By 
inspection of the self-focusing equations, we see that this 
condition is only met if the second term on the right-hand- 
Side is larger than the first. The second term is due to 
the nonlinearity in the dielectric and is thus the focusing 
term while the first term is the diffraction term and if it 
dominates the focusing term , then the beam will tend to 
diverge (i.e. OF 77 A steady state can occur if the two 


terms balance exactiy. 
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4. SELF-FOCUSING IN A MAGNETOPLASMA 


In this chapter, we investigate solutions to the 
equations derived in Chapter 3 which described self-focusing 
i.e. equations {3.12}, {3.14} and {3-15}. These equations 
do not admit an easy analytic solution except under special 
conditions, however a good deal of information can be 


Obtained without resorting to a numerical solution. 


4.1 STEADY-STATE SELF-FOCUSING - SELF-TRAPPED SOLUTIONS 


In this section we investigate self-trapped (or soliton 
like) solutions where the beam propagates in a steady-state 
Manner {i.e. with no convergence or divergence). This 
occurs when the seif-focusing process is exactly balanced by 
diffraction. The ponderomotive and weak heating cases wiil 
be considered first where the nonlinear part of the 
dielectric re is a function of E2, followed by a treatment 


of the strong heating case . 


In the steady-state case we take a— » a, (a, wili 
denote the equilibrium bean size), f(z) =1 and all 
derivatives vanish ({ 4f= We a 2, f=0) ogw ErOMe {3.42} wer obtain 


for the wavenumber shift 
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and from {3.14} we obtain for the beam size 


L 2a 
{4.2} Gat=1mG peeks. 


2 Spor 
Equation {4.1} is the dispersion relation for the nonlinear 
medium for equilibrium propagation. This can be seen more 
clearly if we substitute the expressions for 772 and a, into 


{4.1}, 


Ae 
{4.3} obk ste Bohs o60) - 2G; 4p 


where y (9) is the on-axis density ratio. This dispersion 
relation resembles that obtained for plasma filled 
waveguides where the last term is the correction due to the 
fact that the beam is bounded (i.e. self-trapped), and the 
second term shows the effect of the plasma depletion in the 
presence of the beams. Thus the beam appears to propagate in 


a self made waveguide, or light pipe. 


Equation {4.2} gives the beam size, (as a function of 
the on-axis intensity E?) that is required for self-trapped 
propagation. | With a=agy the beam propagates uniformly 
through the plasma with no change in size. For a given 
aii this gives the critical power required for self- 


trapped propagation (Rx ares) 
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From {4.2} we can obtain a condition for the equilibriun 


beam size, a 


fp! to reach a minimum by setting daf/dE?=0. 


This gives 


{4.5} Egeh =a — ioae/ 5a 


as the condition for the minimum equilibrium beam size a 

to occur. It will be shown later that the existence of 
“1 pet is a result of the limitation on the beam size due to 
diffraction effects and saturation of the nonlinearity. 


‘We now investigate these equilibrium quantities for the 


cases discussed in Chapter 2. 


a) Ponderomotive nonlinearity in a strong magnetic field. 


we treat the case where fB,<<1 so that some analytical 


progress can be made. From {2.19} we have that 


{4.6} ; o(E.) = top. [> Ar Pe 


EG Vion 


once 


so that C3 is parabolic in E. From {4.2} we see that for 


this case 
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The equilibrium beam size for the field free (B=0) case has 


recently been obtained (Max (1976)) and was found to be, 


{4.8} at = acterp|~,Fs 7 | 
lop, , EY 

Thus the equilibrium beam size is much larger for the 
magnetized plasma than for an unmagnetized plasma when 
Ona S 1. The magnetic field case reguires more power to 
achieve self-trapped propagation than the field free case 
because the magnetic field reduces the plasma depletion in 
the beam and hence reduces the ‘'strength' of the _ self- 
focusing effect. A larger beam size is therefore required 


to reduce the diffraction effect (which is proportional to 


1fa2) to balance the weaker selif-focusing. 
The critical power is found from {4.9} to be 


z 


3 
{4.9} heeo= 6620 Mahcmeeme 
Wp. B.%e 


which gives, for CO, laser, B~0.1, bp?hy * 0.1 and T~ 


100eV a Bee = 2x1011 Watts. 


b) Weak heating in a magnetic field. 


From {2.438} we have that 
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which for OK ES<<1, gives the parabolic form 


{4.11} E> we ax E 
Lage 


The equilibrium beam size is 


| L 
aN Ve eG E/ 


{4.12} a 
lop. ae A 


which for OBS << te gives the critical power 


3 r 
{4.13} Tee tave ad Xl Oecs 
hy %, 


For co, laser, Tow) 100eV, by*/is Zn 0-1, this gives a Pe 
HO MW. Thus the ponderomotive force in a strong magnetic 
field requires considerably more power than the heating 
nonlinearity in order to achieve the self-trapped state. 
This is due to the ponderomotive force being inherently weak 
until high intensities are achieved. We thus conclude that 
the heating of the plasma and subsequent density depletion 
is by far the more dominant mechanism for self-focusing for 


laser intensities of interest in present day laser - plasma 
p 


experiments. 


As the heating nonlinearity is not parabolic for 
KEZ~ 1, we can compute eda, from {4.10} using {4.5}... The 
x 


condition for minimum equilibrium beam size to occur is 


u 
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Substitution of {4.14} into {4.12} yields 


{4.15} OE lee a AN 


which is of the order of a few vacuum wavelengths across 
(for & 2g W?)- Similar results have been. found for the 
ponderomotive and weak heating cases in the field-free 
(B, =0) situation. Max (1976) and Sodha et al (1974) find 
Belay C/p, as the typical scale length for these situations 
even though the mechanisms for the nonlinearity are 
different. Also, Kaw et al (1973) find that scale lengths 
of the order of CAsp across the beam are the smallest scale 
Vengths for which the frlimentatton™instability will oecur. 
Max (1976) has suggested that the self-trapped states, which 
are equilibrium states, represent the final stage of the 
filimentation instability because tHe scale lengths are the 


same. We discuss the validity of obtaining such small scale 


lengths across the beam in a later chapter. 


It is interesting at this point to examine the 
dispersion relation {4.3} for a cut-off beam size, ae les 
o 

whereby c2k2=0, and solutions to the wave equation do not 


propagate. This cut-off is 


{4.16} as pd fs 
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occur for beam sizes of the order of a vacuum wavelength. 
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will therefore not exist for Wy ?/s 2<<1. 


Cc) Strong heating in a magnetic field. 


At this point we solve for the equilibrium solutions to 
equations {3.12} and {3.15} for the strong heating case 
Sstudiede sin . 2e3< Setting f(z)=1 and ,£=0 in,y {3315} ,0we 


obtain for the equilibrium radius 


2 r ( 
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wo” oy F'( Xe) 
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The dispersion relation obtained from {3.12} by setting 9,9=0 
is identical to {4.3} except that r! is given by {2.84}. 
From {4.17}, we find that the minimum equilibrium beam size 


Occurs for 


{4.18} %, BS = gf BOVE 


2 


and this gives 


Qa Ci aeteast see Be Saty Rie? ax 
{4.19} plSgye 7 7 bp. 


So again the typical minimum beam size is of the order of 
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C/hw as found in the previous discussion on weak heating in 


P. 
part (b). 


4.2 NONEQUILIBRIUM SELF-FOCUSING - WEAK FOCUSING LIMIT 


We now examine the physical significance of the 
equilibrium beam size by discussing the cases for which the 


function g is parabolic in E for example the 


eo’ 
ponderomotive case for f<<1 or the weak heating case when 
Oo ES <<1. For these parabolic cases, by substituting {4.2} 
into the self-focusing equation {3.14}, the following resuit 


is obtained 
{4.20} SS ( t-4t 
aj 


Thus, if a2>aze then ee and self-focusing will occur 
(i.e. £(Z) will decrease). Self-focusing will not occur if 
a?<a? since 9,570 and the beam will diffract ccntinuously 
Withee ft. ncreasing with) zZ.97So0 a, is the critical beam size 
for which self-focusing will occur and for a given 
intensity, ag determines the critical power, that above 


which, self-focusing Willis ,Occur. Fort ~Bo< Pa es no self- 


focusing will take place. 


We now investigate nonequilibrium solutions for the 
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beam size parameter f£ by introducing the dirtrraction 
distance Ry=ka2, which is a § familiar quantity in linear 
optics. Its significance can be seen by setting the 


rocusing ‘term in {3/14} 2 toezero’ (ive. E2 —p 0) so that it 


becones 
0..3¢2) = halal 
+4 Regs 
a 
with solution 
2 2 
{4.21} Se Ra 


where the. boundary conditions £(0)=1 and af (0) =0 have been 
used. Thus a beam propagating a distance z=R will have 
increased its beam size by a factor of y2’ due to 
dittiraction. We can introduce a similar quatity Re =kaa,, 
which we will cail the ee Wesvomenyernie), distance, and it 
represents the distance the beam has to propagate for the 
width to decrease by a factor of 2, in the absence of 


diffraction. In terms of Ry and R we can write {4.20} as 


BS 


ae 7 ey eras be 
{4.22} ng 762) = = ( R: a ) 


Tros’*self-focusing will, only occur if Es Rae (ve. 1 t 


focusing dominates over diffraction). We recall that {4.22} 
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will only be valid for the cases where ¢ iS. parabolic; gin 


E,- 


The solution of {4.22} with £(0)=1 and 2, £ (0) =0 is 
Sa Pe ee af | ) z 
{4.23} S12 ot Cee. ) e 


SO that. LOT Ry <Rye the beam will converge slowly. 


For the ponderomotive nonlinearity in a strong magnetic 


field, we have from {4.7} 


{4.24} R, = 2kca™ p.<<! 
lup* Po% Ee 


and for the weak heating case 


Re = asc gan ow, BA << | 


FOr, -4a CO laser heated plasma with, 4?/5 0.1, Tr~100 ev 
and 6,~ 0.1, {4.243 and {4.25} say that a self-focusing 
distance of 0.1 cm requires intensities of 101% W/cm2 and 
- 1019 W/cm? respectively, where again we see that the 
ponderomotive nonlinearity requires much higher intensities 
than the weak heating case to give the same effect on 


focusing. 


In Figure 3 we plot {4.23} for various values of 
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Ra/R, - It can be seen how the beam converges more rapidly 
as the self-focusing distance decreases. It must be noted 
that {4.23} is only valid for v3, E2/f2<<1 CEE OE e/a 50 
that as f decreases the approximations will be in jepardy. 
Thus the prediction of catastrophic self-focusing 
(i.e. f —~ 0 as z—> RR / (RE-R?) ) by {4.23} is invalid, and 


we expect new features to appear when f gets smaller. 


4.3 NONEQUILIBRIUM SELF-FOCUSING - OSCILLATORY WAVEGUIDE 


STRUCTURES 


In this section we examine solutions of the self- 
focusing equation for which the non-parabolic form for ry rs 
used. This corresponds to situations where the initial 
intensity is sufficientiy high or where focusing has caused 
the intensity to become very darge such that the 
approximations used in the previous section break down. The 
ponderomotive mechanism is not treated in this section as no 


analytic solutions can be obtained for mw E2>1. 


We begin by considering the weak heating case and 


rewrite equation {3.14} by defining a new function U,(f). 
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where $ is given by {2.58}. Then the self-focusing equation 


{3.14} can be written as 
{4.27} drcadjet- De U(5) 


Equation {4.27} yields, upon integration, the conservation 


law 


{4.28} (2,5)"- GCL eat 


where C is a constant. Using the boundary conditions at z=0 


Eis. | 


{4.29} Ay 
Boo) = £ (@) ATES 


o 


where the last condition is derived from the fact that f3 (9) 
is identified as the initial radius of curvature of the 
incident beam. If the incident beam is a finite plane wave, 
then 9,£(0) =0 = R,—yeo - Thus, evaluating {4.28} at z=0 
with the boundary conditions {4.29} determines the value of 
C as 


{4.30} Ge Se, 
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We can now make use of {4.28} to compute the nmininun 
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value attained by the beam size parameter f. The existence 
of fi, comes from the fact that as the beam focuses, the 
effective intensity of the beam increases, causing the 
second term on the right-hand-side of {3.14} (the focusing 
term) to decrease in magnitude. This decrease in magnitude 
is due to the fact that ¢’ is essentially the derivative of 
the density ratio %, and as the beam focuses, the increased 
intensity causes more and more density depletion in the 
beam; thus ols 0 and ”—P0 as £-> 70. Therefore the 
nonlinear self-focusing saturates at high intensities. At 
the same time as the nonlinearity weakens, the diffractive 
‘torce’) (the RU*f8 term in {3.14}) is increasing. Thus at 
some point the diffraction can dominate the focusing and the 
beam will eventually start to diverge, even though initially 


the beam was converging. 


Using {4.28}, together with the fact that 9, f=0 at 


f=f and taking Rio, we find that E naan satisfies 
{4.31} Wise.) aM) 


We have solved {4.31} for f,;, asa function of KES 3 these 
results appear in Figure 4 and show the saturating behaviour 
of the focusing. Similar results have recently been 


published by Sodha et al (1976). 


The strong heating case (2.3) will now be considered in 


much the same way. We begin by writing the self-focusing 
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equation {3.15} explicitly as 


| } 2 $* ox E 
{4.32} %,5c2)= —— 3) + - we AM 
es Kg eye oe Ce (1426 6,)'" (14% E, ) 


By defining a function Uy (f) as 


uv 
{4.33} VOUS tel Ses ox, £, Ind 
> 2K S* Wy ene yay a Be 
ka” (14.2%, F.)” (14%, &) 


then {4.32} becomes 

{4.34} Og DGB aS) Op ly, F)) 

as we had before. Thus to compute i -, we have to solve 
{4.35} Uz (Smee )2 U20)) 


where we have used the boundary conditions O,f (0) =0 and 
£(0)=1. Thes= solution ~0£84 {4035} 4las8¢a VEunction® of OE, 
appears in Figure 5 where again we see the saturation 


behaviour. 
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5. SOME PHYSICS AND SOME APPLICATIONS OF SELF-FOCUSING 


In this chapter we discuss the previous analytical 
results in terms of some simple physical arguements and draw 
some conclusions based on a physical understanding of self- 
focusing. We also discuss some possible applications of 


these ideas to the laser - solenoid fusion concept. 


5-1 PHYSICAL DISCUSION OF SELF-FOCUSING 


‘The self-focusing process comes down to a simple 
competition between two effects; the tendency of the beam to 
focus due to the density depletion in the beam caused by 
radiation pressure or heating, and the tendency of the beam 
to defocus due to diffraction. The dominant of the two 
effects determines the nature of the beam propagation. A 
balance can also occur between the two effects, as was seen 
in 4.1, where the diffraction of the beam was exactly 
compensated for by the focusing action (of the plasma 
resulting in the propagation of a uniform beam. For a given 
intensity, this balance condition allowed us to determine 
the beam size and hence the power for which this self- 
trapped or equilibriun propagation state would exist. In 
4.2 we saw that for a given intensity, focusing would only 
occur if the beam size 'a' was larger than the equilibriun 


value a This is because a larger beam size reduces the 
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diffractive: "force and: the, focusing, action. is.-able ‘to 


dominate, causing the beam to converge. 


According to the weak focusing solutions of 4.2, once 
the focusing began it was irreversible and catastrophic, 
i.e. the beam would shrink to zero ata point focus, 
causing the intensity to become infinite, which is an 
unphysical, if not uncomfortable conclusion. However, as 
the intensity increases, the validity of the approximations 
used in 4.2 is destroyed and a more rigorous discussion must 
be used. In 4.3 we found that the beam size reached a non- 
zero minimum vaiue, which we can expect for two reasons. 
First, it is well known in optics that the focal spot size 
ofdany opticalcsystem is elimitediethby diffraction mands this 
will prevent the beam from focusing to a point. Second, the 
nonlinearity Saturates at high intensity causing the 
strength of the focusing term to decrease to zero. This can 
be seen in the sketch (Figure 6) of the nonlinear term r) of 
the dielectric. As the focusing depends on the derivative 
of ry With respect to the electric field, then the 
saturation of the nonlinearity at high intensities causes 
$'—30 thus weakening the focusing term. At the same time 
the diffraction effect is getting stronger since the beam is 
getting smaller. Thus, even though the focusing effect may 
indeed dominate initially, the subsequent convergence of the 
beam causes the noniinearity to eventually weaken (saturate) 
and diffraction to increase to the point where it dominates, 


turning the focusing around and thereby spreading the bean. 
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However, as the beam spreads, the diffraction effect begins 
to weaken and we can eventually return to the initial 
Situation where the focusing will again dominate over 
diffraction. Thus an oscillatory behaviour ensues comprised 
of alternate focusing and defocusing. This is depicted-in 


Figure 7. 


The ettects of diffraction and isaturation “of ‘the 
nonlinearity can also explain the existence of the mininun 
equilibrium beam size or es calculated in 4.1. The typical 
size for Seles was ~ Chip for all mechanisms. This size is 
(for by gw) of the order of the free space wavelength of the 
radiation and is a general feature of saturable nonlinear 
nedia (Akmanov et ale. 1966))% This represents the 
diffraction limited beam size for if tg Ree guttiraction ois 
so strong that the nonlinearity cannot balance diffraction 
so that an eguilibrium condition cannot be attained. Even 
if the intensity is increased the nonlinearity weakens 


(saturates) and is less effective in compensating for 


diffraction. 


Dec EOE FECTS OF ABSORPTION ON SELF-FOCUSING 


In the previous analysis of self-focusing (Chapter 4) 
we had assumed that absorption of the laser could be ignored 
which meant that the absorption length, ets waS much 


greater than lengths over which the intensity varied due to 
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focusing (or defocusing). For example if | eas Ray then we 
would expect the effect of absorption to be significant. As 
the absorption length varies as T,* /nz (for? “classical 
inverse-bremsstrahlung) we should solve for self-focusing 
together with absorption in a self-consistent manner i.e. 

any change in intensity brought about by absorption will 
change the plasma parameters which in turn affect the 
absorption. However, this is extremely di ftacuLt 
analytically so we will satisfy ourselves by examining a 


simplified picture and look for gross features. 


We take the absorption coefficient to be a constant and 


replace the intensity E2 b E2e-%? where k is constant. 
iy [<) of ° a 


The weak focusing limit, 4.2 is treated in this way so that 


the self-focusing equation {4.22} becomes 


' Beets 
{5.1} > foe h(a - CC — 
ze S. Ra oe 


The soluticn to {5.1} with no absorption (i.e. Ke =0)@ is 
2 r Noy act 
{5.2} $= ts (re Rie) ee, 
' 5 


from which we can define a focal length Zs (i.e. when f=0) 


or, 
mere 
{5.3} Be = Rak, 
Slnr aio ak 


Now Zp is an order of magnitude estimate for the scale 
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length of axial variation of the beam so that the neglect of 


absorption requires the condition 
7 
{5.4} 2 << Lab’ 


fer {o.4} holds, then for distances’ much “smaller than) (the 
absorption length, e- Kin and {5.1} becomes the equation 
previously solved in 4.2 where absorption was neglected. 
However, if the absorption length is very short, then the 
Magnitude of the self-focusing term in {5.1} decreases 
rapidly with z and if the absorption length is sufficiently 


Short such that 
Lr 
r 
{5.5} Zs a7. a ie 


then for distances larger than a few absorption lengths, 


{5.1} becomes 
{5.6} Ong TCR) ant Ale 
with solution 
{5.7} | S*= Lat t/Rr> 


Thus if the absorption is sufficiently strong, it can cause 
the beam propagation to become diffraction dominated and the 


beam will thus spread, even though initially it may be 
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focusing. Thus, absorption is a defocusing mechanism which 
is to be expected as it weakens the nonlinearity by 
decreasing the laser intensity. Sodha et al (1976) have 
recently solved {5.1} numerically and have demonstrated the 


defocusing nature of absorption. 


In Figure 8 we have plotted the weak (equation {5.2}) 
and strong (equation {5.9}) absorption solutions of {5.1}. 
For intermediate values of the absorption length, solutions 
to {5.1} must lie in the region bounded by the weak and 
strong absorption solutions as is clearly seen in the 
humerical results of Sodha et al . We also note that even 
if the absorption is weak, all beams will eventually defocus 
once the propagation distance has become comparable to the 
absorption length so that self-focusing can only be 


maintained for distances of the order of i Ta 


5.3 SCALE LENGTHS OF SELF-FOCUSING 


In this section we turn our attention to the size of 
scale lengths over which quatities of interest vary, namely 
the beam size and focusing distance (or oscillation 
wavelength). The latter quantity refers to the axial scale 
length of amplitude variation brought about by the focusing 
and defocusing of the beam while the former is a radial 
scale length. In particular we wish to examine the periodic 


propagation structure found in the discussion of the weak 
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and strong heating cases. As the temperature and density 
are functions of the electric field, then the alternate 
focusing and defocusing of the beam in the plasma will cause 
the temperature and density to vary periodically over the 
Same scale length (which we cail WE”) creating an axial 
nonuniform plasma. We wish to estimate this scale length to 


see what axial variation to expect. 


We recall the results of the weak focusing solutions of 


section 4.2 where we could define a focal distance Zs to be 


2 ru 
{5.8} en ce 
Ry as Re 


which for fe Rae could be quite long. However, we are 
interested in solutions for higher intensities (i.e. og EP 
ope oO, E, OF 20: (1))% For these cases we can write the self- 


focusing equations {3.14} and {3.15} as 
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" 2 
R(t BZ 5%) 
j “Fo /2 
5 
{5.10} R, (5) = 
( 14 2%5E,) (14% E, ) 
S* ox; Ey 


where § 7=WpPa2/c?. If we take Rs (£) ian, {5.10} .to be, .given 
by Rs (1) (i.e. f£=1) and treat this as a constant, then 


{5.9} has the simple solution obtained in {4.1} of 


2S 
{5-11} cots [i- Ra Jz, 
Ra) Ry 


and setting f=0 yields, in analogy with {5.8}, a ‘focal’ 


: / 
distance Zz of 


us 


X 
{5.12} zie Ra Kyl) Kt) 
ite ul at) 
Equation {5.12} gives an estimate , albeit a crude one, for 
the scale length of axial variation of the electric field. 
The oscillation length would be 2 os) FOL “EP Or KE, of 
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re ru 
{5.13} 2 ajk gee! te “tk ee 


a ip. 


Although the derivation of {5.13} is not a very exact one, 
it reveals an important property of the focal or oscillation 
Pendgunwine cidat “Chis 2engtne vs of therorder 7ot the initia 
beam size ‘a' for the higher intensities. This can be seen 
more rigorously if, following Max (1976) we consider a beam 
propagating in the equilibrium self-trapped state and 
perturb the solution slightly. This is done by writing the 
right-hand side of the self-focusing equation as a function 
mn {a2,f), and then expanding ¥ about the equilibrium in 


order to linearize the equation. 


We treat the weak heating case here, as the strong 
heating case gives very similar results. Following the 


above procedure, equation {3.14} becomes 


4 
{5.14} ya 562) = AOS F¢2>) 
where 


2 2 
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Ka $? a (1406, E452) ° 
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Then using {5.163}, {5.14} can be written as 


Pred A AS SE] bitc( Sal heiee Cehi ane 2 egenat 
he Kat (V4855) k a, 


which has the solution that satisfies f=1 and 2, £=0 at z=0, 


Lee 
(5.18) 5-|=2(a¢-a") |- cos! za: ( As an Jz] 


eaten pave a, ES 
kata? paige 


Thus the perturbed solutions oscillate with os, GtVen by 


SS pt Ve 
{5.19} Asse = V2 7 ka, ( UBL A) 


As the initial beam size wassa,, then we see again that the 
oscillation scale length is determined, in part, by the beam 
size. ‘This feature also appears in treatments of bean- 
trapping where the density profile is assumed to be given, 
instead of solving for it self-consistently (see Mani et al 
(1974), Feit and Fleck (1976)), where it is found that Ais 
determined by the width of the density profile. As the bean 
size determines the width of the density profile (in the 
steady state) then this is the reason for the connection 


between beam size and A. 
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for beam sizes of the order of the radiation wavelength r\ ’ 
then A~o(A), which violates the assumption that the axial 
Variation of the amplitude is slowly varying (see discussion 
after {3.5}) which helped to simplify the wave equation. 
Thus for the previous discussion to be valid, we require the 
initial beam size '‘a' to be much larger than the radiation 
wavelength, i.e. a>>k-1. Therefore, beams that are of the 


order of alee C/wp. tend to violate our slowly varying 


Pon po 
assumption, so that for ge se our solutions should have 
A, o>k he even for high intensities. This is in contrast 
with the conclusion of Max (1976) whereby she claims that 


Az ko for high intensities, regardless of initial beam 


size. 


5-4 SMOOTHING OF INHOMOGENITIES 


Mani et al (1974) pointed out that the axial variations 
of fenperacure and density that accompany the periodic self- 
focusing wili be smoothed out by axial thermal conduction. 
The conduction aiong magnetic field lines is due primarily 
to electrons and lengths over which conduction keeps’ the 
temperature uniform is given by 
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where @ is the time scale involved. For a1 keV plasma 
with Hee eeOee Cha, hpiew Tener for’“CO1i0r nsec. As the 
heating time scales of interest in fusion reactors are 
nearing Ipsec, then the temperatures and hence density 
should be fairly uniform over distances of tens of 
centimeters. In constrast, rose is oftehestordéerAnottsan few 
beam sizes, and typically a~0.1 - 1.0 cm, so we see that 
Nese <<ih at very high temperatures and thus the 
inhomogeneity in density and temperature dqetsatced By the 


non-uniform intensity should be smoothed out. 


In the radial direction, perpendicular to the magnetic 
field, the density and temperature variations are determined 
by the beam width, and as the beam focuses, the radial width 
of the density and temperature profiles also decrease. 
However, radial thermal conduction will limit the size of 
these profiles. Heat conduction across field lines is 
primarily due to ions because of their larger Larmor radii, 
and the length over which radial heat conduction will keep 


the temperature and density uniform is given by 


{5.21} 4. = jo" Sen eae: 


For Bw 100 kG, LS RORe cm—~3, I,~1 keV then {5.21} gives 
L~ 10h and 103” for &~ ins and lms Cespectively. So for 
CO, radiation (k~!1~ 104), even for short time scales, it is 
unlikely that beam sizes of a~w~k-1 would self-focus as a 


density profile that narrow would not be maintained 
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sufficiently long before conduction smoothed it out. Thus 
for a longer time scale, the perpendicular heat conduction 
can cause the temperature and density profiles to be fairly 
uniform across the extent of the beam, which will tend to 
weaken the focusing effect, causing the periodic focusing to 
Snooth out. In fact, if the’ radial -thermal conduction i is 
very ‘efficient', then this could untrap the beam by not 
permitting a density profile to be maintained. TRUS! 0 
order to cut down on radial thermal conduction, a strong 
magnetic field inside the plasma will probably have to be 
used in order to trap a laser beam over large distances, as 


required for laser - solenoid fusion schemes. 


5-5 APLICATIONS TO LASER SOLENOID FUSION 


We have seen in our discussion of the ponderomotive and 
heating mechanisms that the latter is by far the more 
dominant self-focusing process for intensities of interest 
in laser - solenoid fusion (1011 - 1012 W/cm?2), even in 
strong magnetic fields. However, the focusing of the laser 
beam due to heating and subsequent expansion of the plasma 
causes the intensity to increase by factors as high as 102 - 
103 so that the ponderomotive force may become important in 
these high intensity regions either by causing further 
plasma depletion or by coupling the transverse e.n. 
radiation to longitudinal plasma modes (i.e. stimulated 


scattering processes). Nevertheless we have seen that the 
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heating of the plasma can create a density profile that is 
'favourable' to affect focusing even though the plasma was 


initially uniforn. 


In the foregoing analysis, we have computed the density 
profiles self-consistently, but these are only valid when a 
steady-state condition has been reached. This means that 
these profiles are established on time scales of the order 
of the characteristic heating or phonon transit times (Gor 
a! see introduction). Thus we have not solved the self- 
focusing problem for time scales shorter than these, leaving 
the question of the focusing (or defocusing) of the initial 
portion of the beam unanswered. We can, however glean some 
qualitative information on the development of the initial 
portion of a laser puise by considering a slowly increasing 
pulse. If the rate of increase of the pulse is sufficiently 
slow so that we can apply our weak heating theory as an 
approximation over time scales At SCs » then from the 
self-focusing equation {3.14}, we see that for early times 
in the. pulse, where the intensity is very small, the 
focusing term can be neglected and diffraction initially 
dominates. For later times the diffraction should become 
compensated for by the focusing action as_ the pulse 
intensity begins to increase and at some point the self- 
trapped state should be attained. As the pulse intensity 
increases further, the periodic focusing and defocusing 
discussed in 3.3 will ensue. If we recall that absorption 


is a defocusing mechanism and that for early times the 
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plasma will be relatively cold implying short absorption 
lengths, then this process will tend to cause the early 
portions of the pulse to defocus as well as restrict any 
focusing that may take place, to regions of length less than 


hye These gualitative features appear in Figure 9. 


The qualitative conclusions, as shown in Figure 9, have 
been born out in a recent time dependent numerical treatment 
of the plasma radial dynamics and axial beam propagation by 
Feit and Fleck (1976), where they have seen features similar 


to Figure 9, even for short pulses. 


From the above discussion, we can conclude that the 
initial portion of the laser pulse will diffract out of the 
plasma as sufficient time has not elapsed to establish a 
density profile sufficient to affect beam-trapping. fThus, 
unless there exists a pre-formed, favourable density 
profile, the beam front wiil tend to defocus. This will 
probably also be the tendency of the beam just behind the 
bleaching wave front (the sharp boundary between optically 
thick and thin regions of the plasma column, see Steinhauer 
and Ahlstrom (1975)). However, the beam may produce a 
density minimum ahead of the bleaching front because in the 
regions far behind the front the plasma is being heated and 
is expanding, thereby blowing the magnetic field imbedded ina 
the plasma out as well. This will cause field lines further 
ahead to become distorted and take some plasma with then, 


thereby creating a density minimum ahead as well. Now the 
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disturbance of the magnetic field will propagate at the 
Alfvén velocity Vy=B/T4xnn; so that if the bleaching wave 
propagates at or slower than VU, then the heating and 
subsequent expansion of the plasma (with the magnetic field) 
will tend to create a favourable density profile ahead of 
the front and may tend to focus the beam (or at least weaken 
the diffraction tendency) (Burnett and Offenberger (1976)). 
However, this requires a closer look at the coupled radial 


and axial dynamics of the laser - plasma interaction. 


By taking diffraction into account, we have seen in 
3.3, that even though catastrophic self-focusing does not 
occur, the intensities attained in the focused portions of 
the beam can be very high indeed. This is quite undesirable 
as this puts the intensity into a regime where parametric 
instabilities such as stimulated Brillouin scattering (Drake 
et al (1974)) may be important. Also this strong focusing 
can create nonuniformities in plasma density and temperature 
as discussed earlier. To reduce the strong focusing, we 
could use very large B-fields so that a nearly constant 
radial density profile is maintained thereby weakening the 
focusing tendency. As we saw that the strong focusing will 
not take place until some heating has caused a density well 
gs form, the application of strong B-fields would not be 
reguired until the plasma was ‘hot'. However, near the bean 
front, as the plasma temperature is cooler, high magnetic 
fields would not allow an appreciable density profile to 


form so that beam trapping may not even occur. 
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Thus to control the focusing in some fashion, the 
Magnetic field would have to be increased in the high 
temperature, heated region of the plasma column and not in 
the bieaching front region. So a pre-determined rate of 
increase of B along the axis could help to reduce the strong 
focusing in the hot regions by reducing the amount of plasma 
depletion in the beam. High magnetic fields would also be 
neccessary to prevent diffractive losses at the ‘hot! end of 
the solenoid (the end where the laser is incident) as radial 
heat conduction can smooth out a favourable density profile 
when high temperatures are reached unless a high magnetic 
field is present to reduce this effect. But the use of high 
Magnetic fields is undesireabie from the point of view of 
obtaining high beta plasmas for greater efficiency in fusion 
reactors, so that this is one problem which still has to be 


resolved. 


To avoid the use of high magnetic fields to control 
strong focusing tendencies, the laser energy could be 
distributed over many different laser lines so that the bean 
would be composed of radiation spread over many wavelengths. 
As the focusing distance Zs OL oscillation wavelength Nose 
depend on wavelength (see equations {5.13} Por e(5:319})t> “then 
the use of a multimode laser would tend to smear the 
focusing out since some of the modes would be focusing while 
at the same time others would be defocusing (Offenberger 


(1977)). As opposed to the case where all the laser energy 
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is in a single monochromatic line, giving rise to high 
intensities, any particular mode of a multimode laser would 
have a somewhat lower intensity so that strong focusing 
would tend to be less severe in terms of the high 
intensities that it would produce. This is the result of 
only some of the laser energy being focused into a_e small 
region at any particular point. Also the smearing out of 
the focusing would tend to smooth out inhomogeneities in 


temperature and density caused by nonuniform heating. 


We end this chapter by mentioning those areas related 
to the propagation of laser beams in long plasma columns 
that require further investigation. Srrst, Of vali; ~the 
Stability of the laser beam - plasma system requires 
investigation, Steinhauer (1976) has shown that the beam is 
stable to long wavelength, axial perturbations (long 
compared with the focusing length), however Feit and Maiden 
(1976) claim that short wavelength perturbations can cause 
exponential growth of the beam size. Both of these 


investigations are quite preliminary. 


Also, scattering due to nonlinear processes such as 
stimulated Brillouin or Raman scattering (SBS or SRS) may 
become important if the beam focuses to high intensities. 
Sodha et al (1976) have shown that self-focusing in an 
unmagnetized, collisionless plasma can enhance SRS 


considerably. 


Perhaps the most important piece of work that is 
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required, apart from an actual experiment, is a completely 
self-consistent solution of the coupled, time dependent, 
plasma - optical problem where both radial and axial plasma 
dynamics are taken into account. Also the propagation of 


multimode liasers in plasmas requires investigation. 
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6. CONCLUSION 


In this thesis we have shown that the nature of the 
laser beam propagation depends on the competition between 
the refractive (focusing) effect of the plasma and the 
diffractive (defocusing) tendency of the beam. Asa _ result 
of this competition, the beam will either defract, propagate 
uniformly, or alternately focus and defocus, with the focal 
size limited by diffraction. This behaviour can be produced 
by either a collisionless ponderomotive mechanism or by the 
heating of the plasma. We have studied the effect of an 
externally produced magnetic field inside the plasma and 
Shown that more laser energy is required to produce the same 
effect on the focusing of the beam as compared to the field 
free case, as the plasma depletion is reduced by the B- 
field. It was also shown that regardless of the self- 
focusing mechanism, the smallest equilibrium beam size is of 
the same order, whether or not a magnetic field is present. 
Thus we conclude that it is diffraction, which is common to 
all cases considered, that limits the beam _ size. The 
heating mechanism was found to be a far more dominant 
process than the ponderomotive nonlinearity for Aces powers 
ai interest to solenoid fusion, although the ponderomotive 


force may become important in the focal regions of the bean. 


We have speculated that strong absorption of the laser 


may tend to defocus the beam front for early pulse times 
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while strong radial conduction may flatten the density 
profile at late times (when the plasma is hot) and _ hence 
destroy the focusing mechanisn. A pre-formed density 
Minimum may be necessary to heip focus the initial bean 
front while to counter the effects of the strong focusing 
that can occur later, we have suggested that strong magnetic 
fieids might have to be used or that the laser energy be 
distributed over many laser lines (i.e. use of a multimode 
laser). However, high magnetic fields will probably be 
required at late times to reduce the strong radial thermal 


conduction which could destroy the focusing profile. 
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Table 1. Restrictions for Weak Heating Model 


° Ee<< B, >> 
(eV) (statvolt/cnm) 2 (kG) 


10 1.3x108 7.2 


102 1.3x109 23 
103 - 1443x1910 72 


1.3x1011 230 


(Note: the restriction on the magnetic field for a given 
temperature can be expressed as /=BEn, 1, /BE<<7-7x10Tt Ee, 
which gives #<<1 for n,=7x10!7 cm-3.) 
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Fig. Be Beam size parameter f(z) - Weak focusing limit. 
See equation {4.23}. 
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FEQ.,, fe “hibernate focusing and diefocusing of Laser bean in 
a plasma. 
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c distance into plasma 


Fig. 9. Qualitative time developement of laser propasztion 
in a plasma. (a), and (b) correspond to early times in the 
pulse where diffraction doninates and spreads the beam 
although the spreading in (b) iS somewhat diminished by 
focusing tendencies. (c) and (d) correspond to later times 
where focusing is balanced by diffraction ((c)) or dominates 
diffraction ((d)), although absorption may defocus the bean 
after propagating a short distance into the plasma. Ani (8) 
represents the late times when focusing completly dominates 
diffraction and the absorption lengths are very long. 
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APPENDIX A DERIVATION OF EQUATION {2.45} 


We derive the equilibrium equation for electrons by 
setting the d.c current to zero (i.e. J9=0) . In an axial 
magnetic field there exists both a radial and 6-component of 
J°so we set each component to zero to obtain the egquilibriun 
equation. The same equation is obtained from both 
components so that we will treat only one, the -component 


here. The §$-component of the zero-order (d.c.) current is 
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Setting J° to zero yields the desired equation, {2.45} 
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APPENDIX B DERIVATION OF EQUATION {2.78} 


From {2.76} we have 


Differentiating twice w.r.t. xr yields 
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Setting r=0 in {B2}, we obtain, using {2.73} and {2.77}, 
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so that for XK E=1, this term is 0.2/5 *, and it is safe to 

neglect the eH) term and we then obtain {2.78} from {B3} 
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Grea 
{B 5} Lids wa 
If the last term in {B3} cannot be neglected, then its 
effect of the temperature profile width 9” would be to widen 
it. We can state that the validity of {B5} requires the use 
of sufficiently strong magnetic fields as oo E,—v% and 
% (0) ~—p 9 when B,—p0. Piece e. vecomes, ©too large: the 
assumption of a Gaussian shaped temperature profile is 


probably in doubt as well. 
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